Quasiparticle dispersion of the t-J and Hubbard models 
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The spectral weight A(p,tj) of the two dimensional t — J and Hubbard models has been cal- 
culated using exact diagonalization and quantum Monte Carlo techniques, at several densities 
1.0 < (n) < 0.5. The photoemission (uj < 0) region contains two dominant distinct features, namely 
a low-energy quasiparticle peak with bandwidth of order J, and a broad valence band peak at en- 
ergies of order t. This behavior persists away from half-filling, as long as the antiferromagnetic 
(AF) correlations are robust. The results give support to theories of the copper oxide materials 
based on the behavior of holes in antiferromagnets, and it also provides theoretical guidance for the 
interpretation of experimental photoemission data for the cuprates. 

PACS numbers: 74.20.-z, 74.20.Mn, 74.25.Dw 
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Angle resolved photoemission (ARPES) techniques ap- 
pHed to the high temperature superconductors have pro- 
duced interesting data that introduces important con- 
straints on theories for the copper oxide planes. Re- 
cently, it has been shown Q that the hole-doped com- 
pounds Bi2Sr2CaCu208, Bi2Sr2Cu08, YBa2Cu307, and 
YBa2Cu408 exhibit universal properties likely induced 
by the behavior of carriers in their common Cu02 planes. 
In particular, it has been reported that the quasipar- 
ticle dispersion has a small bandwidth governed by an 
energy scale of the order of the exchange J of the Heisen- 
berg model ('^ 0.15eV). In addition, in the vicinity of 
momenta Y=(0,7r) and X=(7r,0), the dispersion is 
anomalously flat. These results give support to theo- 
retical ideas based on strongly correlated electrons, since 
(i) it is well-established |^ that at half-filling the spectral 
function of a hole in an antiferromagnet contains a sharp 
quasiparticle peak at the top of the valence band spectra 
with a bandwidth of order J, and (ii) careful studies of 
the fine details of the hole dispersion in one band mod- 
els have revealed the presence of flat regions near the X 
and Y points in momentum space. The existence 

of these two features is a direct consequence of the pres- 
ence of strong correlations and antiferromagnetism in the 
cuprates. 

It is reasonable to assume that the behavior of holes in 
systems with long-range antiferromagnetic order will not 
change qualitatively as the density of holes is increased 
away from half-filling, as long as the antiferromagnetic 
correlation length ^af remains large. Theories based on 
this assumption have been proposed. In particular, 
in Refs. |^,^ it was shown that it is possible to reproduce 
many of the anomalous properties of the cuprates, includ- 
ing the presence of a d-wave superconducting state and 
the existence of an optimal doping, with the economical 
assumption that the sharp quasiparticle peak observed at 
half-filling at the top of the valence band remains robust 
as the electronic density decreases to phenomenologically 



realistic values. This assumption (i.e. approximate rigid- 
ity of the quasiparticle dispersion with doping) received 
support from recent calculations addressing the presence 
of "shadow bands" in the cuprates. The rigid band hy- 
pothesis has also been studied by other authors. [|l^ On 
the experimental side, recent ARPES results by Aebi et 
al. |ll] have shown that features induced by the AF cor- 
relations at half-filling are also present at optimal doping. 
Since the closest structure to the Fermi level in A(p, us) 
is likely to dominate the low temperature properties of 
the model, then it is important to establish theoretically 
whether the quasiparticle peaks observed at half-filling 
survive in the presence of a finite density of holes. 

The purpose of this paper is to report results of an 
extensive analysis of the spectral weight for both the 2D 
t — J and Hubbard models using exact diagonalization 
(ED) and quantum Monte Carlo (QMC) methods, sup- 
plemented by Maximum Entropy (ME)) techniques, and 
carried out at several densities. A(p,aj) is shown to con- 
tain a two-peak structure, with dispersing features near 
the top of the valence band dominated by the scale of 
antiferromagnetism J, while a secondary broad structure 
appears at energies of order t. We discuss the range in 
parameter space where this behavior is to be expected, 
and its influence on the physics of carriers in the cuprates. 
However, note that recent QMC results have reported the 
presence of only one PES peak for the Hubbard model at 
both half-filling jl^ and finite hole density. Q We found 
that the disagreement with our present results is avoided 
once the influence of finite temperature effects is consid- 
ered, and a more sophisticated ME method is used. 

The technical details of the present computational 
study, as well as the Hamiltonians of the Hubbard and 
t — J models, are the standard ones, unless otherwise 
stated. In Fig. la, A{p,uj) is shown for the t — J model 
at half — filling and J/t = 0.4 using the ED technique 
applied to 2D clusters with 16 and 18 sites. The combi- 
nation of these clusters allows enough resolution in mo- 
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mentum space to quantitatively analyze the dispersion 
of the main features in the spectral weight. The present 
results have been obtained using approximately 100 itera- 
tions in the standard continued fraction expansion (CFE) 
method to obtain dynamical properties using the Lanc- 
zos technique. j|] However, the figure shows that only a 
small number of poles are dominant. It is clear that near 
the Fermi energy, (oj = 0), there is a robust peak that 
weakly disperses in the scale of the figure. Remnants of 
this low-energy peak exist at momenta (0,0) and (7r,7r), 
in the latter barely visible to the eye (but its intensity and 
position can be easily studied with the CFE approach 
mentioned above). In Fig. lb, the position of the low- 
energy peak is shown with full dots, with the convention 
that the area of the dot is proportional to the intensity 
of the peak. The bandwidth of this sharp quasiparticlc- 
like peak is ~ 0.8t — 2J, in excellent agreement with our 
expectations based on previous ED and Born approx- 
imation [|j calculations. The flat region near (tt, 0) is 
also visible in the figure. From Fig. la it is clear that 
additional PES spectral weight in A(p, uj) is located at 
higher energies |u;|. As discussed before in the literature, 
the strong correlation effects force the hole quasiparti- 
cle to carry only a fraction of the integrated weight, jT^ ] 
and thus the presence of considerable incoherent intensity 
deep in energy is reasonable. A rough estimation of their 
position is shown in Fig. lb (open squares). Q This fea- 
ture is not relevant for the low temperature behavior of 
the model which is dominated by the quasiparticle peak 
at the top of the valence band. 

Before describing the density dependence of our re- 
sults, let us clarify the importance of finite size effects 
in Fig. la, b, as well as the differences between our re- 
sults and those of previous QMC simulations. |^,|l^ To 
address both issues simultaneously, we have carried out 
an extensive QMC simulation of the 2D Hubbard model. 
The results reported here correspond to U/t = 10 (i.e. 
the strong coupling regime where the model should be- 
have similarly to the t — J model), temperature T = t/4, 
and using 10^ sweeps over the entire lattice to reduce 
the statistical errors. Here we use the "classic" ME tech- 
nique. ||l^ This method gives a closer fit to the Monte 
Carlo data than the variant used in Ref. [p|jT^ and there- 
fore resolves more structure. The analytically-continued 
A(p, uj) obtained at half-filling on an 8 x 8 cluster is shown 
in Fig. 2 at several momenta. The results are both quali- 
tatively and quantitatively similar to those obtained for 
the t — J model in Fig. la, b, and also in good agreement 
with ED studies for the Hubbard model, jl^] Two peaks 
in the PES region are clearly identified for all momenta. 
From their position it can be shown that the bandwidth 
of the peak at the top of the valence band is of order 
J, in excellent agreement with our previous discussion. 
The second broader feature observed in the ED study 
of the t — J model is also present in the QMC simula- 
tion results. Studies at larger U/t couplings in the Hub- 



bard model and in the region 0.2 < J/t < 0.8 of the t — J 
model show basically the same features. Then, here it 
is concluded that the qualitative physics of both mod- 
els is very similar in the strong coupling region, where a 
A(p, Lu < 0) with a double-peak structure is observed, as 
properly assumed in previous analytical studies. Jzplp^ 
Let us now discuss our results away from half-filling. 
In the relevant regime of density for the high-Tc super- 
conductors, namely in the vicinity of "optimal doping" 
(n) w 0.85, the QMC-I-ME method at large U/t produces 
stable results only at temperature T = t/2 which is too 
high to resolve the two peak structure even at half-filling. 
Thus, in this density regime only the ED results are re- 
liable. In Fig. 3a, b, ED data at density (n) « 0.88 are 
shown (two holes in the 16 and 18 sites clusters). In 
this case (^af is approximately two lattice spacings. 
The PES results along the diagonal in momentum space 
present structure very similar to that discussed at half- 
filling. The low energy peak is well-defined at all mo- 
menta, even those above the naive non-interacting Fermi 
momentum located near (7r/2,7r/2), and still it disperses 
with a bandwidth of order J. The large accumulation of 
weight at higher energies \uj\ remains localized at w ~ At. 
Then, to the extend that the one band models reproduce 
the physics of the high-Tc cuprates, it is reasonable to 
expect that PES experiments carried out at half-filling 
and near the optimal doping, should produce dispersive 
features of similar intensity and bandwidth. The clear 
similarity between the experimental bandwidth of the 
Bi2212 PES data, and recent results for the insulating 
Sr2Cu02Cl2 compound, provides more evidence for 
the validity of strongly correlated one band models for 
the cuprates. However, it is important to remark that 
while the concrete prediction of our calculations is that 
the bandwidth of the hole carriers is of order J, the par- 
ticular details of the dispersion may dif fer from com- 
pound to compound. For example, it has been recently 
remarked that to reproduce the data for Sr2Cu02Cl2, 
the addition of a small t'-tcrm to the 2D t — J model is 
necessary. Thus, care must be taken when the fine 
details of different compounds at different dopings are 
compared. 

Now consider the inverse photoemission (IPES) {lu > 
0) intensity in Fig. 3a, b. The observed spectral weight 
in the vicinity of (tt, tt) somewhat resembles the distri- 
bution for a non-interacting Fermi system. In principle, 
this effect does not seem reproduced by a rigid band fill- 
ing of the states at half-filling. However, recently Eder 
and Ohta 120] have shown that if proper quasiparticle 
operators |13| are used in the calculation of the spectral 
weight (i.e. operators dressed by spin fluctuations, in- 
stead of bare electronic operators), then the intensity of 
the IPES region is much reduced and the quality of the 
rigid band description of the t — J model appears more 
clearly. This is an important point not much emphasized 
in the literature, namely that the robustness of the rigid 
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band picture in a given model cannot be tested by ana- 
lyzing the removal of "bare" electrons (sudden approxi- 
mation) as produced by a PES experiment, but instead 
"dressed" carriers must be used. Thus, PES and trans- 
port experiments may differ in their predictions if holes 
are heavily renormalized as in the cuprates. 

Fig. 4a, b shows ED results for A(p,a;) using the same 
clusters and coupling as at half-filling, but now reduc- 
ing further the density to (n) « 0.75 and 0.50 (i.e. 4 
and 8 holes in the 16 and 18 sites clusters). In this case, 
through the spin correlations we observed that ^af is less 
than one lattice spacing and thus the influence of AF fluc- 
tuations should be small at these densities. Indeed the 
two-peak structure discussed before at higher densities is 
now difficult to identify. While the broad valence-band 
feature at w ~ 4t remains, only remnants of the AF- 
induced intensity above the naive Fermi momentum can 
be observed. The IPES signal increased its intensity and 
now A(p, w) resembles the behavior of a non-interacting 
^^cospx + cospy" band. 

An interesting detail of Figs, la, 3a and 4a, b, is that the 
intensity of PES weight at p = (tt, tt) changes appreciably 
as the density is varied. This is to be expected since p = 
(tt, tt) is the momentum the most sensitive to the presence 
of AF correlations. In particular, when ^af 0, we 
expect that the PES weight at p = (tt, tt) will be mostly 
transferred to the IPES regime. The presence of PES 
weight at (tt, tt) and (n) = 1 is a direct consequence of 
the AF correlations, and for a paramagnetic background 
A((7r, 7r),cj < 0) should be negligible. 

Summarizing, in this paper the quasiparticle disper- 
sion of the 2D t — J and Hubbard models was analyzed 
as a function of the electronic density. At half-filling, 
A(p, Lu < 0) has a sharp quasiparticle- like peak at the top 
of the valence band with a bandwidth of order J. This 
structure is the relevant one for the low temperature be- 
havior of the models. A second broad feature deeper in 
energy was also identified. As the electronic density de- 
creases, the "two peak" structure remains clearly visible 
as long as the antiferromagnctic correlation length ^af 
is robust. When the AF fluctuations become negligible 
then a crossover exists into a dispersion for the quasipar- 
ticles which resembles a weakly interacting system. For 
realistic values of the coupling, namely U/t = 10, this 
crossover from an antiferromagnctic metal to a paramag- 
netic ground state occurs between (n) = 0.88 and 0.75. 
Then, in the interesting regime for the copper oxide ma- 
terials the AF correlations govern the behavior of the 
spectral weight. The present results give strong support 
to theories of the cuprates based on the behavior of carri- 
ers in an antiferromagnet, and provides information 
about the crossover from a half-filled to a doped system 
that can guide the analysis of ARPES data. 
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Figure Captions 

1. (a) Spectral weight A(p,a;) of the 2D t — J model 
at J/t = 0.4 using clusters of 16 and 18 sites along 



the diagonal in momentum space. The i5-functions 
have been given a width e = 0.25i in the plots; (b) 
position of the two dominant peaks in A(p, lu) as 
a function of momentum. The area of the circles 
is proportional to the intensity of the quasiparti- 
cle peak they represent. The error bars denote the 
width of the peak as observed in Fig. la (sometimes 
to a given broad peak several poles contribute ap- 
preciably). The full squares at lu ^ —At represent 
the center of the broad valence band weight, and 
the area of the squares is not proportional to their 
intensity. 

2. Spectral weight A{p,uj) of the 2D Hubbard model 
obtained with the QMC method supplemented by 
Maximum-Entropy, on an 8 x 8 cluster, U/t = 10, 
and T = t/4. 

3. Same as Fig.l but for density (n) « 0.88 (i.e. two 
holes on the 16 and 18 sites clusters). In (a) the 
PES intensity is shown with a sohd line, while the 
IPES intensity is given by a dotted line. The chem- 
ical potential is located at ut — 0. In (b) the full 
and open circles represent the PES and IPES in- 
tensities, respectively, of the peaks the closest to 
the Fermi energy. Their area is proportional to the 
intensity. 

4. (a) Same as Fig. 3a but for density (n) « 0.75 i.e. 4 
holes on the 16 and 18 sites clusters; (b) Same as 
Fig. 3a but for density (n) w 0.50 i.e. 8 holes on the 
16 and 18 sites clusters. 
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